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Uracil DNA glycosylase (UNG) loss enhances DNA 
double strand break formation in human cancer cells 
exposed to pemetrexed 

LD Weeks 1 , GE Zentner 2 ' 5 , PC Scacheri 2 ' 3 and SL Gerson*' 1 ' 3 ' 4 

Misincorporation of genomic uracil and formation of DNA double strand breaks (DSBs) are known consequences of exposure to 
TS inhibitors such as pemetrexed. Uracil DNA glycosylase (UNG) catalyzes the excision of uracil from DNA and initiates DNA 
base excision repair (BER). To better define the relationship between UNG activity and pemetrexed anticancer activity, we have 
investigated DNA damage, DSB formation, DSB repair capacity, and replication fork stability in UNG +/+ and UNG " 7 " cells. 
We report that despite identical growth rates and DSB repair capacities, UNG ~'~ cells accumulated significantly greater uracil 
and DSBs compared with UNG +/+ cells when exposed to pemetrexed. ChlP-seq analysis of y-H2AX enrichment confirmed fewer 
DSBs in UNG +/+ cells. Furthermore, DSBs in UNG +/+ and UNG" 7 " cells occur at distinct genomic loci, supporting 
differential mechanisms of DSB formation in UNG-competent and UNG-deficient cells. UNG " 7 " cells also showed increased 
evidence of replication fork instability (PCNA dispersal) when exposed to pemetrexed. Thymidine co-treatment rescues S-phase 
arrest in both UNG +7+ and UNG" 7 " cells treated with IC 50 -level pemetrexed. However, following pemetrexed exposure, 
UNG 7 " but not UNG +7+ cells are refractory to thymidine rescue, suggesting that deficient uracil excision rather than dTTP 
depletion is the barrier to cell cycle progression in UNG ~ 7 ~ cells. Based on these findings we propose that pemetrexed-induced 
uracil misincorporation is genotoxic, contributing to replication fork instability, DSB formation and ultimately cell death. 
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Pemetrexed is a multi-target antifolate that inhibits several 
folate-dependent enzymes, leading to a global reduction in 
nucleotide synthesis and cell death. 1 In addition to reduced 
nucleotide biosynthesis, inhibition of thymidylate synthase 
(TS) - the major target for pemetrexed activity - directly limits 
the methyl reduction of deoxyuridine monophosphate (dUMP) 
to form deoxythymidine monophosphate (dTMP) 2-4 The 
resulting imbalance in deoxyuridine and deoxythimidine 
nucleotide pools favors the utilization of dUTP in DNA 
replication and the accumulation of uracil in DNA. 2 

Uracil DNA glycosylase (UNG) is the major glycosylase for 
the removal of uracil from DNA and the initiation of base 
excision repair (BER). 2,5-9 As a mechanism of TS-inhibitor- 
induced cell death (or so called thymine-less death), it has 
been hypothesized that futile cycles of uracil excision by UNG, 
BER pathway repair activity, and dUTP re-insertion contribute 
to the fragmentation of genomic DNA. 4 ' 10 Though uracil 
accumulation and UNG excision activity have been observed 
in cells treated with TS-inhibitors, 11 little evidence supports 



futile cycles of BER as effectors of thymine-less death. If TS- 
inhibitor-induced genomic instability were reliant upon UNG 
excision of uracil, one would anticipate increased TS-inhibitor 
sensitivity when UNG expression is elevated. However, 
pemetrexed-resistant subtypes of human cancer have statis- 
tically higher UNG expression than pemetrexed sensitive 
subtypes. 12 Additionally, overexpression of UNG does not 
sensitize cells to raltitrexed or plevitrexed, a similar agent. 13 
Recently, we have reported that UNG is predictive of 
pemetrexed response in experimental human lung cancer 
models and that loss of UNG expression in human cancer 
sensitizes cells to pemetrexed. 12 These data imply that 
pemetrexed-induced genomic instability is mediated by 
accumulation rather than removal of uracil. The objective of 
the present study was to further investigate the pemetrexed 
hypersensitivity and DNA double strand break (DSB) forma- 
tion observed in UNG _/ ~ cells in order to better understand 
the consequences of deficient uracil excision during peme- 
trexed exposure. 



department of Pathology, Case Western Reserve University School of Medicine, 10900 Euclid Avenue, Cleveland, OH 44106, USA; department of Genetics and 
Genome Sciences, Case Western Reserve University School of Medicine, 10900 Euclid Avenue, Cleveland, OH 44106, USA; 3 Case Comprehensive Cancer Center, 
Case Western Reserve University, 2103 Cornell Road, Cleveland, OH 44106, USA and department of Medicine, Division of Hematology/Oncology, Case Western 
Reserve University School of Medicine, 2103 Cornell Road, Cleveland, OH 44106, USA 

"Corresponding author: SL Gerson, Department of Medicine, Case Comprehensive Cancer Center, 111000 Euclid Avenue, Wearn 151, Cleveland, OH 44106, USA. 
Tel: +216 844 8562; Fax: +216 844 4975; E-mail: slg5@case.edu 

5 Current address: Basic Sciences Division, Fred Hutchinson Cancer Research Center, 1100 Fairview Avenue North, Seattle, WA 98109, USA 
Keywords: base excision repair; uracil DNA glycosylase; double strand breaks; pemetrexed; uracil; thymineless death 

Abbreviations: BER, base excision repair; ChlP-Seq, chromatin immunoprecipitation high-throughput sequencing; CldU, chlorodeoxyuridine; DNA, deoxyribonucleic 
acid; DSB, double strand break; dTMP, deoxythymidine monophosphate; dUMP, deoxyuridine monophosphate; IdU, iododeoxyuridine; TS, thymidylate synthase; 
TSS, transcription start sites; UNG, uracil DNA glycosylase 
Received 02.9.13; accepted 10.10.13; Edited by H-U Simon 



UNG loss increases pemetrexed-mediated DSBs 
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Herein, we show that despite equivalent rates of prolifera- 
tion, DLD1 human colon cancer cells with stable knockout of 
UNG (UNG _/ ~) are hypersensitive to pemetrexed-induced 
DNA DSB formation and cell death. Using y-H2AX chromatin 
immunoprecipitation followed by high-throughput sequencing 
(ChlP-Seq) we observed significantly greater number of y- 
H2AX-associated genomic loci (DSB) in pemetrexed UNG ~'~ 
cells compared with UNG +/+ cells. Additionally, ChlP-Seq 
data reveal distinct patterns of DSB susceptible loci in UNG +/+ 
and UNG _/ ~ cells, suggesting disparate mechanisms of DSB 
formation in UNG-proficient and -deficient cells. Using the 
complimentary approaches of CldU/ldU DNA labeling and 
PCNA dispersal experiments, we show evidence of replication 
fork instability in pemetrexed-treated cells. UNG +/+ and 
UNG _/ ~ have similar basal level DSB repair protein 
expression and DSB repair capacities, yet UNG _/ " cells 
show delayed recovery from pemetrexed-induced DNA 
damage, including DSBs. Based on these findings we propose 
that the enhanced pemetrexed cytotoxicity in UNG _/ " 



cells is indicative of the genotoxicity of misincorporated uracil, 
which persists at sufficient quantities to stall DNA replication 
triggering cell death. 

Results 

Loss of UNG expression hyper-sensitizes human cancer 
cells to pemetrexed. Previously, we showed that UNG 
expression predicts sensitivity to pemetrexed in experimental 
models of human lung cancer. 12 As we were interested in the 
impact of UNG expression on pemetrexed response in 
human cancer cells, we evaluated sensitivity to pemetrexed 
inUNG +/+ and UNG " /_ colon cancer cells. These UNG ~'~ 
cells were engineered as described 14 and lack expression of 
both nuclear and mitochondrial UNG isoforms (Figure 1a). 
Loss of expression of UNG abolished in vitro cutting activity 
in UNG/APE oligonucleotide cutting assay (Figure 1b), 
suggesting that in the absence of UNG, uracil is not efficiently 
removed from DNA. 





Pemetrexed Concentration (nM) 



Figure 1 Loss of UNG enhances pemetrexed sensitivity in DLD1 human colon cancer cells, (a) Western blot of UNG nuclear (top band, 39 kDa) and mitochondrial (bottom 
band, ~ 36 kDa) in UNG + ; + and UNG ~ ; ~ DLD1 human colon cancer cells, (b) UNG cutting activity assay using either purified enzyme (APE and UNG + APE lanes, 1 U of 
each) or 2.5 ^g of whole-cell extract from UNG +/+ and UNG _/ ~ cells, (c) Proliferation rates were determined by cell counting with trypan blue exclusion of non-viable cells 
over a period of 0-96 h using a hemocytometer. Colony survival assays for UNG +/+ and UNG _/ ~ cells treated with pemetrexed (d), raltitrexed (e) and cisplatin (f) for 
10-days. UNG _/ " cells are 10-fold more sensitive to pemetrexed, P< 0.0001 and 4.3-fold more sensitive to raltitrexed, P< 0.0001. Sensitivity to cisplatin was not 
significantly different for the two cell lines, P = 0.1 7. UNG FLAG cells are UNG ~ 1 ~ cells that were transfected with a FLAG-tagged UNG vector to induce re-expression of UNG. 
UNG FLAG cells showed resistance to pemetrexed and raltitrexed similar to UNG +/+ cells (g) Western blot of apoptotic markers cleaved PARP and cleaved caspase 3 in 
UNG +/+ and UNG _/ ~ cells treated for 24 h with 0-100 nM pemetrexed. (h) Percentage of Annexin V-positive pemetrexed-treated UNG +/+ and UNG _/ ~ cells after 24-h 
exposure to 0-100 nM pemetrexed. (i) Tumor volumes of xenograft tumors grown in NODSCID mice that were treated for 5 consecutive days with pemetrexed (150 mg/kg, 
intraperitoneal injection) 
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UNG expression is coordinated with DNA replication. 7 ' 8 ' 15-17 
This relationship prompted us to evaluate the impact of 
UNG loss on cellular proliferation rates for UNG +/+ and 
UNG" 7 " cells (Figure 1c). An identical number of cells 
(1 x 10 5 ) were seeded in 60-mm culture dishes and incubated 
for 24-96 h. Following incubation, cells were stained with 
trypan blue and the total number of viable cells per sample 
was determined using a hemocytometer. This simple experi- 
ment yielded identical proliferation rates for UNG +/+ and 
UNG" 7 " cells (P= 0.292; Figure 1c). Clonogenic survival 
assays revealed that UNG" 7 " cells were 10-fold more 
sensitive to pemetrexed (IC 50 UNG +/+ = 205.8 nM and IC 50 
UNG" 7 " =22.0nM, P<0.0001, Figure 1d). UNG" 7 " cells 
displayed cross-sensitivity to the TS inhibitor, raltitrexed 
(P< 0.0001; Figure 1e), but not the cross-linking agent, 
cisplatin (P=0.17; Figure 1f). The impact of pemetrexed on 
cell killing was further analyzed by cleaved PARP measure- 
ments using western blot (Figure 1g) and Annexin V 
measurements using flow cytometry (Figure 1h) at 24-hours 
in cell treated with 0-100nM pemetrexed. Results indicated 
that UNG" 7 " cells were significantly more sensitive to 
pemetrexed-induced apoptosis compared with UNG +7+ 
cells. Re-expression of UNG (UNG FLAG ) rescued pemetrexed 
and raltitrexed sensitivity in UNG" 7 " cells (Figures 1d-f), 
confirming that TS-inhibitor sensitivity in UNG " 7 " cells is due 
to UNG loss. 

To evaluate the impact of UNG loss on tumor cell growth 
in vivo, UNG +7+ and UNG" 7 " tumors were xenografted 
subcutaneously into NODSCID mice (Figure 1i). Treatment 
with five daily consecutive intraperitoneal (IP) injections of 



pemetrexed (150mg/kg) resulted in tumor quadrupling time 
of 7.83 ±0.88 days in UNG +7+ tumors compared with 
12.54 ±0.42 days in UNG" 7 " tumors, P<0.001. At day 20 
pemetrexed-treated UNG-/- tumors were 61% smaller 
than untreated UNG-/- tumors while UNG +7+ tumors 
were only 26% smaller with pemetrexed treatment compared 
with untreated controls (Figure 1i). Together these data 
indicate a significant increase in the antitumor effect of 
pemetrexed in UNG-deficient tumor cells in vitro and in vivo. 

Increased DNA DSB formation in UNG 1 cells during 
pemetrexed exposure. The formation of DNA DSBs is a 
long observed consequence of TS inhibitor exposure and 
thymine-less death, 10 ' 18,19 yet the mechanism of DSB 
formation has remained ill defined. It has been hypothesized 
that futile cycles of uracil excision, BER and dUTP reinsertion 
contribute to DNA fragmentation and DSB formation in cells 
treated with TS-inhibitors. 4 ' 10 However, this hypothesis lacks 
experimental support 13 and fails to explain the pemetrexed 
hypersensitivity observed in UNG" 7 " cells (Figure 1). 
We therefore evaluated DNA DSB formation in our isogenic 
UNG +7+ and UNG" 7 " cells to better understand the role of 
UNG in mediating DSB formation and thymine-less death 
during pemetrexed exposure. Western blots for the DNA 
DSB surrogate marker, y-H2AX, show significant induction in 
UNG" 7 " cells compared with UNG +7+ cells (Figure 2a). 
Importantly, even when the two cell lines are treated 
with equally cytotoxic concentrations of pemetrexed 
(IC 50 , UNG" 7 " 25 nM and UNG +7+ 200 nM), induction of 
7-H2AX protein expression appears greater in UNG" 7 " cells 
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Figure 2 Increased DNA DSB formation in UNG _/ ~ cells treated with pemetrexed. (a) Western blot for y-H2AX was performed using whole-cell extracts from UNG +/+ 
and UNG _/ ~ cells treated with 25 nM pemetrexed for 0-72 h to compare sensitivity of the two cell lines to pemetrexed-induced DSB formation, (b) Western blot for y-H2AX 
was performed using whole-cell extracts from UNG +/+ and UNG _/ ~ cells treated with IC50-level pemetrexed alone or in combination with 10^M supplemental thymidine 
for 24 h. UNG +/+ cells were treated with 200 nM pemetrexed and UNG _/ ~ cells were treated with 25 nM pemetrexed. (c and d) Single-cell gel electrophoresis (comet) 
assay was performed under neutral conditions for UNG +/+ and UNG _/ " cells treated with IC50-level pemetrexed alone or in combination with supplemental thymidine. 
At least 50 comet tail lengths were measured using NIH Image J software from two independent experiments for a total of > 100 cells per drug treatment. Representative 
images are shown in (c) and the mean and S.E.M. for this experiment is plotted in (d). Pemetrexed treatment caused 1.85-fold increase in comet tail length in UNG _/ " 
compared with UNG +/+ cells, and P<0.05. Thymidine rescued sensitivity to pemetrexed in UNG +/+ (*P<0.001) and UNG _/ " (**P<0.004) cells 
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(Figure 2b, 'P' lanes). As a complimentary experiment, 
we measured DNA DSBs using neutral single-cell gel 
electrophoreses (comet assay). Pemetrexed treatment 
increased comet tail lengths in UNG _/_ cells by 1.85-fold 
compared with UNG +/+ cells, P<0.05 (Figures 2c and d). 
By providing a salvage pathway for dTTP synthesis, 
supplemental thymidine should limit the effects of TS 
inhibition. Thymidine supplementation rescued IC 50 peme- 
trexed-induced y-H2AX induction (Figure 2b) and comet tail 
formation (Figure 2c) in both DLD1 UNG _/ " and UNG +/+ 
cells, suggesting that in both cell lines dUTP/dTTP nucleotide 
pool ratio drives DNA DSB formation. Of note, the increased 
formation of DNA DSBs in UNG~ /_ cells occurs despite 
equivalent expression of DNA DSB repair proteins and 
DSB repair capacities in UNG +/+ and UNG _/ " cells 
(Figure 3). 

To further compare the mechanism of DSB formation and 
thus cell death inUNG +/+ andUNG -7- cells, we performed 
7-H2AX ChlP-seq to compare the genomic distributions of 



DSBs in pemetrexed-treated UNG +/+ and UNG" 7 " cells. 
Pemetrexed treatment resulted in a 5.75-fold increase in 
distinct y-H2AX-bound genomic regions in UNG" 7 " cells 
compared with UNG +7+ cells. Moreover, of the y-H2AX- 
enriched sequences identified in UNG" 7 " cells, only 
0.0025% of them overlapped with y-H2AX bound sites in 
pemetrexed-treated UNG +7+ cells (Figure 4a). Aggregate 
plots of 7-H2AX peaks confirm that enrichment is significantly 
elevated in treated compared with untreated samples and 
reaffirm a more robust enrichment of y-H2AX in UNG " 7 " cells 
(Figure 4d). When a detailed analysis of genomic features of 
y-H2AX-enriched sequences was performed, DSBs were 
observed to map more frequently to transcription start sites 
(TSS), CpG islands, DNAse-hypersenstive clusters and 
putative origins of replication in pemetrexed-treated UNG" 7 
" cells compared with UNG +7+ cells (Table 1 ). In all, our data 
indicate that UNG" 7 " cells accumulate significantly greater 
number of DSBs in response to pemetrexed and differential 
chromosomal distribution of y-H2AX enrichment in the 
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Figure 3 Comparison of DNA DSB repair capacity in UNG +/+ and UNG _/ ~ cells, (a) Western blots comparing basal expression of proteins involved in DNA DSB repair 
in UNG +/+ and UNG _/ " cells, (b) Schematic for HR repair assay. Briefly, transient transfection of ISCE1 linearized pDR-GFP was used to measure homologous 
recombination efficiency. pDR-GFP contains two tandem but inactive GFP repeats, one of which contains an ISCE restriction enzyme cutting site. GFP fluorescence is 
activated by homology-mediated deletion between the two repeats during functional HR to generate a functional GFP. (c) Schematic for NHEJ repair assay. Briefly, pE-GFP- 
Pem1-Ad2 has a non-functional GFP interrupted by Pern 1 intron into which there is inserted an adenovirus exon sequenced flanked by ISCE1 restriction enzyme cleavage sites. 
ISCE1 plasmid linearization generates non-complimentary ended DSB and removes the AD2 sequence. NHEJ allows juxtaposition of the two GFP fragments to reconstitute 
fluorescence, (d) Representative flow cytometry dot blots from untreated UNG +/+ and UNG _/ ~ cells transfected with an empty vector, ISCE1 linearized pDR-GFP or 
ISCE1 linearized pE-GFP-Pem1-Ad2. pdSRed was co-transfected in each instance as a transfection efficiency control. Similar percentages of GFP + /dSRed + cells were 
observed for UNG +/+ and UNG _/ ~ cells for each transfection suggesting similar capacity for BER 
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Figure 4 y-H2AX ChlP-seq in UNG +/+ and UNG _/ - cells. Cells were treated with IC 50 -level pemetrexed (UNG +/+ , 200 nM and UNG _/ -, 25 nM) or with 20 /M 
cisplatin for 24 h and subsequently processed for ChlP-seq as described in the Materials and Methods section. Cisplatin was utilized as a control DNA-damaging agent that 
forms DNA double-strand breaks, (a-c) Venn diagrams of the ChlP-seq peaks generated comparing signal overlap for pemetrexed treated UNG +/+ and UNG _/ ~ cells (a); 
pemetrexed treated and cisplatin treated UNG +/+ cells (b) and pemetrexed-treated and cisplatin-treated UNG _/ ~ cells (c). (d) Midpoint coordinates for each y-H2AX peak 
were used to obtained signals ± 1 kb for the coordinate (within a 25-bp window). The average signal in each window for each peak shown for UNG +/+ pemetrexed treated 
(top left); UNG +/+ cisplatin-treated (top right); UNG _/ " pemetrexed-treated (bottom left) and UNG _/ " cisplatin-treated cells (bottom right) 



genome suggests disparate mechanism of DSB formation in 
UNG-proficient and -deficient cells. 

Delayed recovery from S-phase arrest in pemetrexed- 
treated UNG" 7 " cells. Pemetrexed-treated UNG" 7 " cells 
accumulate 40-fold more genomic uracil compared with 
UNG +/+ cells. 20 Additionally, there are data to suggest that 
heavily uracilated DNA alters DNA-protein interactions due 
to perturbations in basal genomic methylation, which may 
negatively impact replication and transcription. 21,22 Putative 
origins of replication in human cancer cells were identified 
as common sites pooled from published nascent DNA 
sequencing of K562 and MCF7 cells. 23,24 In ChlP-seq 
experiments, we observed more significant association of 
7-H2AX enrichment with putative origins of replication 
(Table 1). We therefore surmised a mechanistic link 
between UNG loss, DNA DSB formation and replication 
instability. 

Both UNG +/+ and UNG " 7 " cells accumulate in S-phase in 
response to 24-h exposure to IC 50 -level pemetrexed 
(Figure 5a). Phosphorylation of chkl (ser345) and cdc2 



Table 1 Characteristics of y-H2AX-enriched sequences 
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(tyr15) also suggested S-phase arrest (Figure 5b). Interestingly, 
UNG" 7 " but not UNG +7+ cells show delayed recovery 
from S-phase arrest when pemetrexed is withdrawn, 
compared with UNG +7+ cells (Figure 5a). Expression of 
replication-activated checkpoint proteins and arrest in S 
phase with delayed recovery suggest pemetrexed-induced 
replication instability. Supplemental thymidine co-treatment 
blocked S-phase accumulation in pemetrexed-treated UNG +7+ 
and UNG" 7 " cells (Figure 5c), suggesting S-phase arrest is 
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Figure 5 Delayed recovery from cell cycle arrest in pemetrexed-treated UNG cells, (a) Cell cycle progression was monitored using PI staining and flow cytometry. 
UNG +/+ and UNG _/ ~ cells treated with IC 50 level pemetrexed (UNG +/+ , 200 nM; UNG _/ ~ 25 nM) both show evidence of S-phase accumulation after 25 h of exposure. 
When allowed to recover in drug-free media, UNG +/+ cells are capable of resuming cell cycle (return of G1 peak in + 48 h sample) while UNG _/ ~ cells are not. (b) Western 
blot from UNG +/+ and UNG _/ ~ cells treated for 0-72 h with 25 nM pemetrexed showing increased sensitivity in UNG _/ ~ cells to pemetrexed-induced phosphorylation of 
chkl and cdc2, S-phase checkpoint kinases, (c) Cell cycle histograms for cells treated with supplemental thymidine. Cells were treated with 1 0 thymidine for 24 h (No Pern, 
+ Thy), co-treated with IC 50 level pemetrexed and thymidine for 24 h (24 h, + Thy) or treated with IC 50 level pemetrexed for 24 h and allowed to recover for 48 h in media 
containing 10 /M thymidine (48 h post, +Thy) 



related to dUTP/dTTP pool imbalance. However, culturing 
cells in medium containing supplemental thymidine after 
pemetrexed exposure accelerates recovery from S-phase 
arrest in UNG +/+ cells but not UNG _/ " cells (Figure 5d). The 
inability of thymidine to restore cell cycle progression in 
UNG _/ ~ cells after pemetrexed is withdrawn suggests that 
DNA damage caused by misincorporated dUTP is irreversible 
in the context of deficient uracil removal. 

The association of H2AX phosphorylation, S-phase arrest 
and replication fork collapse has previously been studied in 
cells treated with ionizing radiation, nucleoside analogs and 
antimetabolites, 25-27 suggesting that H2AX is phosphorylated 
during inhibition of DNA synthesis. We therefore evaluated 
the stability of DNA replication fork in UNG _/ " cells. 
Intracellular replication fork progression is typically studied 
using nucleotide incorporation assays (Figure 6a). Using this 
approach in a flow cytometry-based pulse chase assay, we 
observed a decrease in post-pemetrexed treatment IdU 
incorporation in UNG" 7- cells (Figure 6b) While this is 
consistent with decreased replication fork progression in 
pemetrexed-treated UNG _/ ~ cells, these data alone could 
simply indicate that misincorporated uracil impedes IdU 
incorporation because they compete for nucleotide insertion 
sites. Therefore, as a complimentary experiment, we assessed 
the dispersal of the replication processivity factor, PCNA, from 
chromatin using both flow cytometry (Figure 6c) and western 
blot (Figure 6d). Similar experiments have been used to 
analyze collapsing replication forks in etoposide and hydro- 
xyurea-treated cells. 28 ' 29 In these assays, decreased PCNA 
staining of S-phase-gated cells or chromatin is indicative of 



replication fork instability or fork collapse. Our data show 
increased PCNA dispersal in pemetrexed-treated UNG _/ ~ 
cells compared with UNG +/+ cells, suggesting that 
pemetrexed causes replication fork instability in UNG-deficient 
cells (Figures 6c and d). 

Discussion 

We have previously shown that loss of UNG expression 
sensitizes human cancer cells to pemetrexed and that UNG 
expression predicts pemetrexed sensitivity in experimental 
models of human cancer. 12,20 Here, we analyzed the DNA 
damage response to pemetrexed in UNG +/+ and UNG _/ ~ 
DLD1 human colon cancer cells to better understand the role 
of UNG in the mechanism of pemetrexed-induced DNA DSB 
formation and cell death. We observed that loss of UNG 
hyper-sensitized DLD1 cells to pemetrexed in vitro and in vivo 
despite equivalent proliferation rates. This hypersensitivity is 
associated with increased replication fork instability and DNA 
DSB formation, despite an equivalent capacity for DNA DSB 
repair in the two cells. 

The formation of DNA DSBs in cells treated with TS 
inhibitors has been studied for decades, 10 ' 30-35 yet the 
precise role of uracil misincorporation and UNG excision of 
uracil in the mechanism of DNA DSB formation and cell death 
is not clear. The dominant futile cycle hypothesis proposes 
that DSBs arise as a result of continuous cycles of uracil 
excision, BER and uracil re-insertion. Experimental evidence 
from other labs 13 ' 32 ' 36 and the data presented herein indicate 
that futile cycles of BER do not adequately explain thymine- 
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Figure 6 Intracellular replication fork instability in UNG ~ 1 ~ cells, (a) Schematic for CldU and IdU pulse labeling of pemetrexed-treated cells. Briefly, cells were incubated 
for 45 min with CldU (50 /M) prior to treatment with IC 50 -level pemetrexed (UNG +/+ , 200 nM, UNG _/ " 25 nM) for 24 h. Cells were then allowed to recover in drug-free 
media for 0-48 h prior to incubation for 45 min with IdU (50 fM). Cells were then stained with fluorescent antibodies that recognize CldU (rat-antiBrdU-FITC) and IdU (mouse- 
antiBrdU-APC). Representative data from flow cytometry sorting of CldU/ldU-labeled cells is shown in (b) where the xaxis is CldU (FITC) and the y axis is IdU (APC). (c) Cells 
were treated with IC 50 -level pemetrexed for 24 and 48 h and subsequently stained with PI (xaxis) to and FITC-labeled PCNA antibody (y axis), (d) Cells were treated as in 
(c) and subsequently incubated in 1% formaldehyde for crosslinking and chromatin extraction for western blots of chromatin-bound PCNA and histone H3 (loading control) 



less death. Overexpression of UNG, which should exacerbate 
futile cycles of UNG, does not enhance TS-inhibitor sensitivity. 13 
UNG loss does not cause compensatory upregulation of the 
other DNA glycosylases capable of uracil excision, 12 so it is 
unlikely that futile cycles of BER are initiated in UNG _/_ cells 
treated with pemetrexed. The relatively fewer DNA DSBs 
observed in UNG +/+ cells compared with UNG _/ ~ cells 
treated with equally toxic concentrations of pemetrexed 
implies that, at least in the models examined, UNG activity 
limits rather than promotes pemetrexed-mediated DNA DSB 
formation and cell death. Based on these data, we bring forth a 
novel hypothesis for the mechanism of thymine-less death in 
UNG-deficient cancer cells. In our model, uracil accumulates 
at a critical level near replication origins, stalling DNA 
replication fork progression and leading to fork collapse, 
DNA DSB formation and cell death (Table 2). 

Several levels of experimental evidence support this model. 
First, we have observed substantial accumulation of genomic 
uracil, 20 coordinated with replication fork instability and DNA 
DSB formation in UNG" /_ cells. Otherwise isogenic UNG +/+ 
cells are largely protected from these effects, suggesting 
uracil is intrinsically cytotoxic. The altered methylation status 
and secondary structure of heavily uracilated DNA may 
explain this toxicity. 21 ' 22 ' 37 An alternative explanation for 
uracil-mediated DSB formation in UNG _/_ cells is increased 



endonuclease cleavage in uracil-containing DNA, 38,39 yet the 
activities of endonucleases with this activity have not been 
assessed in pemetrexed-treated cells. 

Secondly, we have noted persistence of S-phase arrest in 
pemetrexed-treated UNG _/ " cells with concomitant activa- 
tion of intra-S phase checkpoint proteins. UNG +/+ cells still 
undergo S-phase arrest and cell death, albeit at significantly 
higher pemetrexed concentrations than UNG _/ ~ cells. While 
genomic uracil does not accumulate 20 in UNG competent 
cells, dTTP levels are still aberrantly low during TS-inhibitor 
exposure. Low dTTP may limit allosteric regulation of 
ribonucleotide reductase, resulting in elevated dATP and 
decreased dGTP nucleotide pools. 40 In UNG +/+ cells, 
S-phase arrest and cell death in response to pemetrexed 
may be a consequence of global nucleotide pool imbalance. 

Thymidine salvage prevents S-phase arrest but has limited 
capacity to reverse it in pemetrexed-treated UNG _/ " cells. 
Correction of dUTP/dTTP ratios is therefore ineffective after 
uracil has been misincorporated in cells with reduced uracil 
excision activity. Literature suggests that replication fork 
reversal is limited during thymine starvation, while fork 
reversal readily occurs during dNTP depletion or ribonucleotide 
reductase inactivation. 41 Lack of replication fork reversal 
during thymine deprivation may explain the failure to recover 
from cell cycle arrest in UNG _/ ~ cells. 
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Table 2 Possible pathways to DSB formation and cell death in pemetrexed-treated cells 



Possible route to cell death 



Problem with hypothesis 



UNG +/+ cells 
Stalled replication due to nucleotide pool aberrations 
caused by reduced dTTP 
Endonuclease cleavage of uracil and/or AP sites 

Futile cycles of BER 



UNG~ 



cells 



Uracil-mediated replication fork stall and collapse 
Uracil-mediated endonuclease cleavage 



Does not explain differential sensitivity in UNG +/+ and UNG 1 cells. 

No endonucleases with uracil activity have been studied in pemetrexed treated cells. 
Products of APE incision are rapidly targeted by activities of Pol/? and downstream BER. 
UNG +/+ cells are more resistant than UNG _/ ~ cells, suggesting excision of uracil 
protects from pemetrexed cytotoxicity. 



Currently there is no available method to directly interrogate various degrees of uracil 
base substitution within living cells. 

No endonucleases with uracil activity have been studied in this context. 



Lastly, distinct patterns of genomic y-H2AX enrichment 
(DSBs) identified in this study suggest disparate mechanisms 
of DSB formation in UNG-competent and -deficient cells. 
In UNG _/ " cells, enrichment of y-H2AX occurred more 
frequently in accessible chromatin; TSS, DNAse-hyper- 
sensitivity regions, and origins of replication tend to be 
nucleosome-depleted and more readily accessible to 
proteins. This feature is consistent with the negative impact 
uracilated DNA has on DNA-protein interactions. 21 ' 22 
Of interest, our analysis indicates that 36.7% of UNG~ /_ 
and 12.67% of UNG +/+ y-H2AX enriched sites overlap with 
putative human origins of replication, 23,24 suggesting that UNG 
excision of uracil protects origins of replication from peme- 
trexed-induced DSB formation. It has been suggested that 
continued replication in thymine-starved cells triggers destabi- 
lization of origins of replication resulting in DSBs at sites of 
origin firing. 32,36 The fact that UNG-deficient cells appear to be 
more sensitive to this mode of DSB formation is consistent with 
prior observations of coordination of UNG-initiated BER and 
the observation of UNG complexed at replication foci. 6-8,17 

ChlP-seq of y-H2AX to localize DSBs has been previously 
performed. 42 It has been suggested that only a small fraction 
of 7-H2AX foci represent actual DSBs 43,44 Given this caveat, 
we have offered comet assay as corroborating evidence of 
DSB formation. Additionally, over-expression of UNG and 
H2AX in cancer and co-expression of these genes in 
pemetrexed-resistant cancers (Supplementary Figures 1 
and 2) 45-47 reflect the association of UNG and H2AX with 
DNA replication. The deposition of endogenous unpho- 
sphorylated H2AX at replication origins, sub-telomeres, and 
TFS coincides with sites of endogenous DNA damage and 
7-H2AX enrichment. 42 We took steps in our experimental 
design to confirm that the y-H2AX enrichment observed was 
not due to altered replication rates or endogenous DNA 
damage. DSBs were not increased in untreated UNG" /_ 
cells, which had identical proliferation rate as UNG +/+ cells. 
Moreover, <1% overlap between UNG +/+ and UNG" 7- 
y-H2AX-enriched sites suggests different mechanisms of 
direct DSB formation in the two cell lines. 

Combined, the observation of fewer DSBs, reduced 
predilection of DSBs for origins of replication, lack of observed 
replication fork instability, and S-phase arrest that is ade- 
quately rescued with supplemental thymidine in UNG +/+ 
cells suggests pemetrexed-induced toxicity is more related to 
dTTP depletion than to dUTP misincorporation. In contrast, 



UNG cells, which accumulate genomic uracil, have 
increased DSBs that localize to replication origins, show a 
significant degree of replication instability and have irrever- 
sible S-phase arrest; pemetrexed toxicity is more related to 
dUTP misincorporation (Table 2). 

Together these findings highlight the importance of genomic 
uracil accumulation in mediating cellular response to peme- 
trexed. Significant enhancement of pemetrexed cell killing in 
UNG" 7 " cells suggests that targeting UNG may have 
therapeutic impact in human cancer. Additionally, our data 
suggest in pemetrexed-treated cancer cells, cell death is 
mediated by direct toxicity of uracil accumulation and by 
nucleotide perturbations induced by dTTP depletion, with 
genotoxic uracil accumulation becoming the dominant mode 
of cell killing as UNG activity declines. Apart from our studies 
using pemetrexed, studies of changes in UNG expression have 
reported little impact on cellular sensitivity to TS-inhibitor- 
induced cell killing. 13,48 In contrast, one study has reported 
diminished DNA replication rates and increased apoptosis in 
UNG-deficient cells even in the absence of drug treatment. 49 
We speculate that variations in sensitivity among previously 
studied UNG deficient cells may be explained by differences 
in induced genomic uracil accumulation in the TS/inhibitor 
models investigated. Indeed, the sensitization index for 
UNG" 7 " cells to raltitrexed was far smaller than the 
sensitization index for UNG" 7 " cell to pemetrexed 
(Figure 1), suggesting that UNG expression and activity has 
a greater role in pemetrexed versus raltitrexed response. 
To address this, we propose future quantitative studies that 
compare genomic uracil and nucleotide pool levels in 
UNG" 7 " and UNG +7+ cells treated with equally toxic 
concentrations of various TS inhibitors. 

Materials and Methods 

Cell lines and reagents. Pemetrexed was purchased from LC Laboratories 
(Woburn, MA, USA). Thymidine, raltitrexed and cisplatin were purchased from 
Sigma Aldrich (St. Louis, MO, USA). Cells were maintained in complete DMEM 
(10% FBS, 2mM L-glutamine) at 37 °C in a 5% C0 2 incubator. DLD1 cells were 
purchased from ATCC and expanded upon delivery into numerous vials of low 
passage cells for cryopreservation. Cell line characterization by ATCC is conducted 
through short tandem repeat (STR) typing. Re-authentication was not performed. 
DLD1 UNG _/ ~ cells 14 were a gift from Sanford Markowitz. 

Western blots. Protein extracts (30 ^g) were resolved by SDS-PAGE, 
transferred to PVDF membrane (Millipore, Billerica, MA, USA) and detected by 
incubations in primary antibodies followed by HRP-linked secondary antibodies at 
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manufacturer's specifications. Chromatin-bound proteins were extracted from 
formaldehyde cross-linked (1%) cells using the Pierce Chromatin Prep module 
(Thermo Pierce, Rockford, IL, USA). Antibody sources: UNG-23936 (39 kDa band, 
nuclear UNG), Nbs1, DNA ligase IV, Ku70, BRCA1, MRE11, and PCNA (Abeam, 
Cambridge, MA, USA); Tubulin (Calbiochem, via Millipore, Billerica, MA, USA); 
7-H2AX, cleaved caspase 3, and Histone-H3 (Millipore); cleaved PARP, XRCC4, and 
secondary antibodies (Cell Signaling, Danvers, MA, USA); and Rad51, p-chk1, chkl, 
p-cdc2, cdc2 (Santa Cruz, Dallas, TX, USA). 

UNG activity assay. UNG activity was measured using 40-mer oligodeox- 
ynucleotide duplex containing a single U:A mispair as described. 12 Products were 
resolved by denaturing polyacrylamide gels (20%), prior to visualization with a 
Typhoon 9200 fluorescence imager (Amersham Bioscience via GE Healthcare, 
Pittsburgh, PA, USA) and UNG activity was calculated (percentage of cut band 
relative to total oligonucleotide). 

Determination of cell proliferation rate/double time. Cell counting 
with trypan blue exclusion determined cellular proliferation rates. Briefly, 1 x 10 5 
cells were seeded in 60-mm cell culture dishes and allowed to grow in typical cell 
culture medium for 24-96 h. At the time of the harvest, cells were stained with 
trypan blue, counted using a hemocytometer and the number of viable cells was 
recorded to generate a growth curve. 

Colony survival. Colony survival was assessed by methylene blue staining of 
colonies formed after 10-day exposure to the drugs indicated. Only colonies 
containing >50 cells were counted. 

Annexin V staining. Cells were stained with Annexin-V FITC and 
propidium iodide (PI) for evaluation of apoptosis by flow cytometry according 
to the manufacturer's protocol (BD Pharmingen, San Diego, CA, USA). At the 
time of harvest, cells were stained with 5 fi\ of Annexin V FITC reagent and 
10/il of PI (5/zg/ml) in 1 x binding buffer (10 mM HEPES, 140 mM NaOH, 
2.5 mM CaCI 2 , pH 7.4) for 30 min in the dark. Apoptotic cells were determined using 
a Coulter flow cytometer (EPICS-XL-MCL). Cell death calculations are the sum of 
early apoptotic (Annexin V positive, Pl-negative) and late (Annexin V-positive, 
Pl-positive) cells. 

Cell cycle and PCNA dispersal assays. PI staining of methanol-fixed 
cells was used for cell cycle determinations. For PCNA dispersal assays, 
FITC-labeled PCNA antibody (PCNA-FITC, Abeam) was added for PCNA 
detection. Both uni-parameter (Pl-only) and dual-parameter (PI + PCNA-FITC) 
analyses were performed on a Beckman Coulter flow cytometer (Brea, CA, USA) 
(EPICS-XL-MCL). Cell cycle histograms (PI) and PCNA dot plots (PCNA-FITC) 
were deconvoluted from > 20 000 events using Flow Jo software (TreeStar, INC, 
Ashland, OR, USA). 

Tumor growth assays in xenograft mice. Tumor cells (5 x 10 6 ) 
in early passage were injected into bilateral flanks of female NODSCID mice 
(~5 weeks old). When tumor volumes reached 100mm 3 , mice were divided into 
control and treatment cohorts (n = 6). Mice bearing tumors were treated with 
pemetrexed (150mg/kg) or 100^1 sterile PBS (control) by daily IP injection for 
5 consecutive days. Tumor measurements were taken every 2 days and the 
response was quantified by calculating the tumor volume. 

Comet assay under neutral conditions. Cells were processed for 
comet tail formation under neutral comet assay conditions according to the 
manufacturer's protocol for cell lysis and single-cell gel electrophoresis (Trevigen, 
Gaithersburg, MD, USA). Tail lengths were recorded for at least 50 comets on two 
separate slides (~100 cells per treatment) using ImageJ software. 

CldU/ldU pulse labeling and flow cytometry. Cells were incubated for 
45 min in media containing 50 ^M chlorodeoxyuridine (CldU) prior to drug 
treatments. Following treatment, cells were incubated for 45 min with iododeox- 
yuridine (IdU, 50 fM). For analysis, fixed cells were incubated in DNA denaturation 
solution (4N HCI, 0.1% Triton X-100 in PBS) for 15 min, followed by neutralization 
washes with 0.1 M Na 2 B 4 0 7 buffer. Cells were then incubated for 30 min at 37 °C 
in 1 mg/ml RNAseA in PBS. Each sample was labeled using rat anti-BrdU-FITC 
(for detecting CldU) and mouse anti-BrdU-APC (for detecting IdU) antibodies and 
counterstained with PI. Uptake of CldU and IdU was determined using a BD LSR II 



flow cytometer. Dot plots of cycling cells (cells in G1/S/G2M) were deconvoluted 
from univariate analysis of PI staining. Bivariate dot plots of CldU/ldU content were 
obtained for cycling cells using Flow Jo software. 

DNA DSB repair assays. DNA DSB assays were performed as described. 50 
Briefly, cells were transiently transfected with either circular pDRGFP or ISCE1- 
linearized pDRGFR for homologous recombination (HR) assays or with circular 
pE_GFP-PEM1-Ad2 or ISCE1 -linearized PE-GFP-PEM1-Ad2 for non-homologous 
end joining (NHEJ) assays. DNA repair plasmids were obtained from Elisa Tichy 
(University of Cincinnati) and plasmid stocks were expanded using bacterial 
transformation. Plasmid pDsRed express (Clontech, Mountainview, CA, USA) was 
co-transfected in each instance as a transfection efficiency marker. Transfection of 
2.5 x 10 5 cells in six-well culture dishes was carried out using Opti-MEM media and 
oligofectamine according to standard laboratory procedures. In all, 5 fig of pDsRed 
Express N1 and 40 ^g of construct DNA were utilized for transfection. Seventy two 
hours post transfection, cells were trypsinized and suspended in cold PBS containing 
0.1% BSA for analysis by flow cytometry. Samples were analyzed using a BD LSRII 
(BD Biosciences, San Jose, CA, USA). Dot plots for DsRed and GFP were acquired 
for > 20 000 events using Flow Jo software. 

y-H2AX ChIP sequencing. ChIP was performed as described, 51 with 10 /*g 
of rabbit anti y-H2AX, p-S139 (Abeam, ab2893) per ChIP. Sequencing libraries were 
prepared as described, 51 with the exception of the library amplification step. Instead 
of Phusion DNA polymerase, Expand High Fidelity Plus (Roche Diagnostics, 
Indianapolis, IN, USA) was used. The following parameters were used for thermal 
cycling: 5 min at 94 °C, 30s at 65 °C and 60s at 72 °C followed by a 7 min final 
elongation at 72 °C. Libraries were sequenced at the Case Western Reserve 
University Genomics Core Facility. The following unique read numbers were 
obtained: WT untreated yH2AX ChIP, 28089742; WT cisplatin-treated yH2AX 
ChIP, 22505752; WT pemetrexed-treated yH2AX ChIP 30494966; UNG _/_ 
untreated yH2AX ChIP 36823368; UNG _/ ~ cisplatin-treated yH2AX ChIP 
38155972; UNG" /_ pemetrexed-treated yH2AX ChIP 52013817. ChlP-seq data 
generated in this study were deposited into the Gene Expression Omnibus under 
accession number GSE50692. After removal of duplicate reads, ChlP-seq data were 
aligned to the hg18 human genome assembly with Bowtie, 52 discarding reads with 
> 1 reportable alignment. Peaks were called using MACS. 53 For MACS analysis, 
the pemetrexed- or cisplatin-treated ChIP was used as the treatment file and the 
untreated ChIP was used as the control file. Fold enrichment cutoffs used were 4 for 
WT cells and 6 for UNG _/ " cells. Peak coordinates are given in Supplementary 
Data 1. 

To determine the overlap of yH2AX with putative replication origins, peaks were 
converted to the hg19 human genome assembly with liftOver (http://genome. 
ucsc.edu/cgi-bin/hgLiftOver). Lists of putative origins of replication determined by 
nascent DNA sequencing in K562 and MCF7 cells 23,24 were combined and regions 
<10bp apart were collapsed into a single putative origin. Lists of peaks were 
intersected with the condensed list of origins with the GFF-overlap feature of the 
ChlP-seq Tool Set (http://havoc.genomecenter.ucdavis.edu/cgi-bin/chipseq.cgi), 
with peaks < 200 bp apart considered overlapping. The UCSC table browser was 
used to overlap yH2AX peaks with various regulatory features. The proportions of 
yH2AX peaks overlapping putative origins and other regulatory features in 
pemetrexed- and cisplatin-treated WT and UNG _/ " cells were compared by % 2 
contingency test. 

Statistical analysis. Data are presented as the mean±S.E.M. of at least 
three independent experiments. Significance, assigned for P-values<0.05, was 
determined by unpaired two-tailed Student's f-test with standard software 
(GraphPad Prism, San Diego, CA, USA). 
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